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The  present  work  describes  the  optimisation  of  a  compact  steam  reformer  for  light  fuel  oil  and  diesel 
fuel.  The  reformer  is  based  upon  a  catalytically  coated  micro  heat  exchanger  that  thermally  couples 
the  reforming  reaction  with  a  catalytic  combustion.  Since  the  reforming  process  is  sensitive  to  reaction 
temperatures  and  internal  flow  patterns,  the  reformer  was  modelled  using  a  commercial  CFD  code  in 
order  to  optimise  its  geometry.  Fluid  flow,  heat  transfer  and  chemical  reactions  were  considered  on 
both  sides  of  the  heat  exchanger.  The  model  was  successfully  validated  with  experimental  data  from 
reformer  tests  with  4kW,  6kW  and  lOkW  thermal  inputs  of  light  fuel  oil.  In  further  simulations  the 
model  was  applied  to  investigate  parallel  flow,  counter  flow  and  cross  flow  conditions  along  with  inlet 
geometry  variations  for  the  reformer.  The  experimental  results  show  that  the  reformer  design  allows 
inlet  temperatures  below  773  K  because  of  its  internal  superheating  capability.  The  simulation  results 
indicate  that  two  parallel  flow  configurations  provide  fast  superheating  and  high  fuel  conversion  rates. 
The  temperature  increase  inside  the  reactor  is  influenced  by  the  inlet  geometry  on  the  combustion  side. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Fuel  cell  systems  comprising  of  PEM  fuel  cells  and  steam  reform¬ 
ers  have  high  potential  as  stationary  or  mobile  auxiliary  power 
units.  In  particular,  the  use  of  liquid  fossil  fuels  with  an  existing 
infrastructure  like  fuel  oil  and  diesel  fuel  promises  an  early  eco¬ 
nomic  breakthrough  of  the  technology. 

The  implementation  of  steam  reformers  in  the  power  range  up 
to  20  kW  for  domestic  fuel  oil  or  diesel  fuel  as  a  catalytically  coated 
plate  heat  exchanger  or  micro  channel  heat  exchanger  has  advan¬ 
tages  over  other  reactor  concepts.  In  comparison  to  a  shell  and 
tube  heat  exchanger  the  micro  channel  and  plate  heat  exchanger 
provides  higher  heat  flow  densities.  A  technical  application  for 
methane  steam  reforming  based  on  a  plate  heat  exchanger  with 
wavy  plates  was  developed  by  Gritsch  [1].  The  reformer  has  a 
volume  of  0.7 1  and  provides  up  to  1 5  kW  hydrogen  at  reactor  tem¬ 
peratures  up  to  T=  1 173  K.  Kolb  et  al.  [2]  developed  a  micro  channel 
heat  exchanger  for  steam  reforming  of  diesel  fuel.  Operating  tem¬ 
perature  is  T=  1 073  K.  The  reactor  with  a  volume  of  0.25 1  generates 
2  kW  hydrogen.  On  the  burner  side  of  the  reactor  anode  exhaust 
gas  is  catalytically  burnt.  Another  application  for  diesel  fuel  steam 
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reforming  is  a  reactor  with  S-shaped  micro  channels  developed  by 
Thormann  [3].  The  reactor  is  designed  for  a  fuel  input  of  0.2  kW. 
Hydrogen  is  supplied  to  the  burner  side. 

The  steam  reformer  examined  in  this  work  is  designed  as  a 
catalytically  coated  micro  channel  heat  exchanger.  Therein,  the 
endothermic  steam  reforming  process  of  light  fuel  oil  is  heated  by 
the  exothermic  catalytic  combustion  process.  The  reactive  flows 
are  conducted  through  parallel  channels.  The  steam  reformer  was 
built  up  as  a  full  scale  reactor  and  a  small  scale  reactor.  The  full  scale 
reactor  with  a  volume  of  0.68 1  and  a  fuel  input  of  1 0  kW  fuel  (hydro¬ 
gen  output  of  12kW)  was  used  for  system  testing.  In  the  present 
study  the  reactor  was  scaled  down  to  1  kW.  The  small  scale  reactor 
is  based  on  the  same  heat  exchanger  is  used  to  perform  a  detailed 
test  program.  Within  another  experimental  investigation  variants 
of  the  scaled  down  reactor  are  applied  to  a  catalyst  screening  and 
geometry  optimisation.  In  the  framework  of  reactor  development 
and  improvement  CFD  computations  of  the  reactor  are  performed 
and  the  results  will  be  discussed  in  detail  in  this  paper. 

2.  Theory  and  modelling 

2.1.  Steam  reforming 

Fuel  oil  and  diesel  fuel  consist  of  a  large  number  of  higher  hydro¬ 
carbons.  A  number  of  investigations  have  shown  that  the  steam 
reforming  of  higher  and  lower  hydrocarbons  can  be  modelled  by 
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Nomenclature 

a,  b ,  c,  d 

exponents  of  reaction  rates 

a* 

coefficient  for  calculation  of  effectiveness  factor  [14] 

C 

molar  concentration  of  species  (molm-3) 

Ea 

activation  energy  (J  mol-1 ) 

ko.fei 

rate  constant  of  reforming  reaction  (mol  kg 
(catalyst)-1  s-1  (molm-3)-(a-b)) 

k0 * 

rate  constant  of  reforming  reaction  corrected 
by  catalyst  activity  (mol  kg-1  (catalyst)  s-1 

(molm-3)-(a-b)) 

ks 

rate  constant  of  combustion  reaction  (mol  kg 
(catalyst)-1  s-1  (Pa)-(c+d)) 

h 

rate  constant  of  backward  shift  reaction  (mol  kg 
(catalyst)-1  s-1  (molm-3)-2) 

kf 

rate  constant  of  forward  shift  reaction  (mol  kg 
(catalyst)-1  s-1  (molm-3)-2) 

Kp 

equilibrium  constant  of  shift  reaction 

P 

partial  pressure  of  species  (Pa) 

R 

universal  gas  constant  (J  mol-1  K-1 ) 

ri,r2,r3 

reaction  rates  (mol  kg  (catalyst)-1  s-1 ) 

reff 

observed  rate  of  reaction  (mol  kg  (catalyst)-1  s-1 ) 

rs 

rate  of  reaction  at  exterior  catalyst  surface  condition 
(mol  kg  (catalyst)-1  s-1 ) 

S,  So 

density  of  active  sites  on  catalyst  of  used  (partially 
deactivated)  and  fresh  catalyst  (kg  (catalyst)-1 ) 

X 

conversion  of  fuel 

Xr 

conversion  of  fuel  at  reference  point  conditions 

X\ ,  X2 ,  X3 ,  X4  conversion  of  fuel  at  variation  point  conditions 

a 

catalyst  activity 

V 

effectiveness  factor 

sulphur  content  of  fuel  (wt.  ppm) 

0* 

normalised  Thiele  modulus  [14] 

a  steam  reforming  reaction  to  carbon  monoxide  and  hydrogen.  In 
literature  two  types  of  equilibrium  reactions  have  been  proposed: 
methanisation  or  methane  steam  reforming  and  CO  shift  reaction 
[4-7].  Since  the  temperatures  in  the  reactor  exceed  973  K  the  for¬ 
mation  of  methane  is  not  supported  by  equilibrium.  According  to 
the  studies  of  Thormann  [3]  the  steam  reforming  process  can  be 
modelled  by  the  following  two  reaction  steps: 

CnHm  +  nH20  nCO  +  (n  +  m/2)  (AH°98  <  0)  (1) 

CO  +  H20  -b-  C02  +  H2  ( AH!?98  =  -41  kj  mol-1 )  (2) 


The  steam  reforming  kinetics  of  heptane  on  a  Ni/MgO  catalyst  was 
examined  by  Tottrup  and  Rostrup-Nielsen  [8,9].  Thormann  deter¬ 
mined  a  power  law  reaction  rate  for  hexadecane  on  a  ruthenium 
based  catalyst  [3  ].  The  corresponding  rate  equation  which  was  used 
in  this  work  for  reaction  (1 )  is  given  by: 


n 


=  k  iC; 


C„H, 


Q-b 

i  h2 


(3) 


/<!  =  k0  exp 


(4) 


The  reaction  rate  constant  k0  is  adjusted  in  subject  to  the  activity 
of  the  catalyst  coating. 

Assuming  that  CO  shift  reaction  is  not  limited  by  kinetic  effects, 
an  equilibrium  reaction  rate  can  be  used  for  reaction  (2): 


r2  =  kfCcoQ^o  ~  kbCH2Cco2 


(5) 


kb  = 


7CP 


(6) 


Kp  =  4290.1  exp  ((-^  -  0.001 )  1.4664)  (7) 

The  equilibrium  constant  I<p  was  fitted  in  the  temperature  range 
between  673  K  and  1273  K. 


2.2.  Combustion 


There  are  two  types  of  reactions  which  can  occur  in  a  catalytic 
combustion  reactor:  homogenous  and  heterogeneous  (catalytic). 
Homogenous  reactions  occur  in  the  bulk  phase  while  heteroge¬ 
neous  reactions  occur  on  the  catalyst  surface  [10].  A  reaction 
scheme  composed  of  four  reaction  steps,  two  heterogeneous  and 
two  homogenous  reactions,  was  purposed  for  a  methane  catalytic 
combustion  reactor  model  by  Hayes  and  Kolaczkowski  [10].  The 
gas  phase  reactions  consist  of  an  oxidation  reaction  of  methane  to 
CO  and  oxidation  of  CO.  The  heterogeneous  reactions  were  defined 
by  the  oxidation  reactions  of  methane  and  CO  to  C02.  Applying  this 
reaction  scheme  to  the  catalytic  combustion  of  higher  hydrocarbon 
leads  to  unsolved  problems  since  reduced  reaction  mechanism  of 
homogenous  combustion  can  not  reproduce  the  point  of  ignition 
for  all  conditions.  Beyond  that  Pfefferle  and  Pfefferle  [11]  pointed 
out  that  the  catalyst  may  also  promote  homogeneous  combustion 
in  the  boundary  layer  by  production  of  gas-phase  free  radicals. 
For  this  reason  the  combustion  is  modelled  by  one  overall  het¬ 
erogeneous  reaction  rate  reproducing  reaction  characteristics  from 
homogenous  to  heterogeneous  combustion: 

CnHm  +  (n  +  m/4)02  nC02  +  (m/2)H20  (A H°98  >  0)  (8) 

Yao  [12]  investigated  the  oxidation  of  alkanes  over  noble  metals. 
The  kinetics  of  the  reaction  was  treated  with  a  simple  empirical  law 
model: 


r3  -  fe3Po2Pc„Hm 


(9) 


The  given  values  for  a  platinum  based  catalyst  and  butane  is  used  in 
the  implemented  model.  The  constant  is  adjusted  by  experiments 
performed  on  a  monolithic  catalyst. 


2.3.  Mass  transfer 

In  the  present  case  the  reaction  rates  at  the  catalyst  surface  are 
influenced  by  film  diffusion  and  diffusion  in  the  porous  catalyst 
with  a  thickness  of  25  p,m.  However,  precalculations  show  that  the 
limitation  of  heat  transfer  in  the  catalyst  does  not  exist.  To  consider 
the  diffusion  in  the  film  by  CFD  computation  a  proper  resolution  of 
the  channels  is  necessary.  The  mass  transfer  in  the  porous  catalyst 
can  be  modelled  by  two  different  methods.  The  first  option  is  to 
determine  the  concentration  profiles  in  the  porous  media  by  for 
example  the  dusty  gas  model  [13]  which  describes  the  mass  trans¬ 
fer  in  the  catalyst.  Due  to  the  low  catalyst  thickness  the  applied 
model  was  derived  from  the  second  method:  the  effectiveness  fac¬ 
tor  concept.  In  this  case  the  catalyst  phase  is  not  considered  and 
the  influence  of  internal  mass  transfer  on  the  reaction  rate  can  be 
described  by  an  intrinsic  effectiveness  factor  ij  which  is  given  by: 


Gottifredi  et  al.  [14]  developed  a  method  for  estimating  the  effec¬ 
tiveness  factor  for  many  kinetic  expression,  geometry  and  site 
activity  distribution.  Applying  this  model  to  the  present  case  the 
intrinsic  effectiveness  factor  rj  is  calculated  using  the  expression 
proposed  by  [14]: 

r\  =  [0*2  +  exp(-a*0*2)]-0'5 


(11) 


M.  Grote  et  al.  /  Journal  of  Power  Sources  196(2011)  9027-9035 


9029 


Catalyst 

surface 

Symmetry 


Reformer 

Plate 

Burner 


Fig.  l.  Reactor  model. 


Fig.  2.  Steam  reformer  with  superheater. 

2.4.  CFD  computation 

The  models  described  in  the  previous  chapters  are  implemented 
in  the  commercial  CFD-code  ANSYS®  CFX  11.0  [15].  In  the  frame¬ 
work  of  reactor  optimisation  numerical  calculations  of  the  reactor 
are  performed. 

By  neglecting  heat  losses  at  top  plates  it  is  possible  to  reduce  the 
reactor  model  to  one  plate  with  two  half  channels  and  symmetry 
plane  (see  Fig.  1). 

The  reactions  are  considered  by  implemented  specific  species 
sources  at  the  catalyst  surfaces.  Surrogate  fuel  tetradecane  was 
applied  to  all  computations  representing  fuel  oil. 

3.  Experimental 

3.1.  Steam  reforming  reactor  1 0  kW 

The  steam  reforming  reactor  was  designed  as  a  catalytically- 
coated  plate  heat  exchanger  to  couple  the  endothermic  SR  reaction 
with  an  exothermic  catalytic  combustion  (see  Fig.  2).  Fig.  3  shows 
the  solid  phase  temperature  measurements  points  installed  in  a 
central  plate.  The  generated  reformate  has  to  be  cooled  down 
behind  the  reformer  in  order  to  provide  optimum  inlet  temper¬ 
atures  for  the  downstream  high  temperature  shift  reactor  (FITS). 


Burner 


Burner 


Therefore  another  heat  exchanger  was  integrated  into  the  com¬ 
ponent  in  order  to  utilise  the  reformate  heat  to  superheat  the 
steam  before  entering  the  reformer.  In  addition  to  the  FITS  further 
gas  clean-up  stages  were  installed  downstream  of  the  reformer  to 
achieve  the  CO-concentrations  required  for  polymer  fuel  cells.  Fuel- 
mixture  generation  on  the  combustion  side  was  carried  out  in  a  cool 
flame  reactor.  Therein  fuel  is  injected  by  a  nozzle  and  mixed  with 
air.  Low  temperature  oxidation  reactions  in  the  temperature  range 
between  570  K  and  720  K  provide  the  required  heat  for  evaporation. 

In  a  500  h  long-term  test  with  light  fuel  oil  (sulphur  con¬ 
tent  £s  =  5  wt.ppm),  the  successful  operation  of  the  fuel  processor 
was  demonstrated.  By  implementing  periodic  catalyst  regeneration 
cycles,  the  initial  performance  of  the  reactor  could  be  achieved  even 
after  500  h  of  operation  generating  a  hydrogen  yield  of  over  72  vol.% 
(dry). 

3.2.  Investigations  on  a  small  scale  reactor 
3.2.1.  Materials  and  methods 

The  small  scale  reformers  (SSR)  are  used  in  order  to  optimise 
the  full  scale  reformer  (FSR).  Similar  to  the  full  scale  reformer  they 
are  designed  as  catalytically  coated  plate  heat  exchangers  coupling 
thermally  the  reforming  reaction  with  a  catalytic  combustion.  Opti¬ 
misation  tests  of  the  SSR  are  carried  out  with  respect  to  the  flow 
field  geometry  and  the  reformer  catalyst.  To  guarantee  a  similar 
behaviour  of  small  scale  and  full  scale  reformers,  the  down  scaling 
is  accomplished  by  reducing  the  number  of  heat  exchanger  plates 
and  maintaining  the  plate  geometry.  For  the  SSR  the  thermal  power 
input  is  Pth  =  1  kW  resulting  in  a  gas  hourly  space  velocity  (GFISV) 
equal  to  25,000  h- 1  at  a  steam  to  carbon  ratio  (SCR)  equal  to  4.  Simi¬ 
lar  to  the  FSR  the  SSR  includes  one  plate  for  solid  phase  temperature 
measurements  (see  Fig.  4).  In  addition,  temperatures  are  measured 
at  the  reformer  inlet  and  outlet  and  burner. 

For  the  tests  a  fuel  mixture  of  Shellsol  A 1 00  (20  wt.%)  and  Shell- 
sol  D  100  (80  wt.%)  was  used.  The  diesel  surrogate  obtained  had  a 
chemical  formula  of  C14.37H27.40  and  is  comprised  of  the  fractions 
summarised  in  Table  1. 

Operating  the  SSR,  water  and  fuel  are  evaporated  separately 
before  being  mixed  at  around  673  K  and  being  superheated  rapidly 
to  a  temperature  of 923  K  in  order  to  prevent  carbon/coke  formation 
upstream  of  the  reformer.  The  reformer  was  heated  by  hydrogen 


Table  1 

Diesel  surrogate  composition. 


Fraction 

Value 

Unit 

Method  of  analysis 

Paraffins 

48 

wt.% 

GC 

Naphthenes 

32 

wt.% 

GC 

Aromatics 

20 

wt.% 

SMS  2728 

Benzenes 

5 

wt.  ppm 

GC 

Sulphur 

<2 

wt.  ppm 

SMS  1897 

9030 
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Fig.  5.  Plate  temperatures  at  set  temperature  1023 1<  and  SCR  equal  to  4. 


combustion  on  the  combustion  side.  To  control  the  operating  tem¬ 
perature  of  the  reformer,  measuring  point  no.6  on  the  temperature 
measurement  plate  (see  Fig.  3)  was  found  to  be  the  point  with  the 
highest  temperature  on  the  centre  line. 

In  the  presented  test  series  the  FSR  flow  field  was  maintained 
for  the  SSR  and  the  reference  catalyst  from  the  FSR  was  applied. 
For  the  reformer  testing  a  test  program  was  executed  varying  the 
SCR  between  3  and  4  at  set  temperatures  1073K,  1023K,  and 
973  K.  Operating  points  were  maintained  for  2  h  each  to  guaran¬ 
tee  steady  state  operation.  The  reformate  gas  was  continuously 
analysed  measuring:  hydrogen  (FI2)  concentrations  by  thermal 
conductivity  detector,  carbon  monoxide  (CO),  carbon  dioxide  (C02 ), 
and  methane  (CFI4)  concentrations  by  infrared  analysis,  and  uncon¬ 
verted  hydrocarbons  (CxFIy)  by  flame  ionisation  detection  (FID). 

4.  Results 

Fig.  5  shows  the  solid  phase  temperatures  measured  in  the 
reformer  operation  at  set  temperature  of  1 023  K.  The  temperatures 
are  decreasing  with  rising  distance  from  the  reformer  inlet  since 
hydrogen  combustion  takes  place  close  to  the  inlet  of  the  reformer 
due  to  its  high  reaction  rate.  High  temperatures  occur  on  the  side 
line  since  an  impinging  stream  is  provoked  by  the  burner  side  inlets. 

Fig.  6  describes  dry  reformate  gas  concentrations  after  conden¬ 
sation  at  278  K  for  reformer  outlet  temperatures  between  923  K 
and  1 073  K.  Flydrogen  yields  of  approximately  70  vol.%  are  achieved 
decreasing  slightly  for  higher  reformer  temperatures.  The  tested 
steam  reforming  catalyst  also  promoted  the  water  gas  shift  reac¬ 
tion  ( WGS).  Evidently,  the  decrease  of  the  hydrogen  yield  goes  along 
with  a  decrease  of  the  carbon  dioxide  and  an  increase  of  carbon 
monoxide,  revealing  that  the  WGS  is  constrained  at  higher  temper¬ 
atures.  The  residual  hydrocarbon  concentration  is  decreasing  with 
high  temperatures,  indicating  that  the  reforming  reaction  only  per- 


Fig.  6.  Dry  reformate  gas  concentrations  at  SCR  equal  to  4. 


Reformer  outlet  temperature  /  K 

Fig.  7.  Residual  hydrocarbons  and  conversion  rate  at  SCR  equal  to  4. 

formed  sufficiently  at  temperatures  from  1013  K  at  the  reformer 
outlet  (see  Fig.  7).  The  fuel  conversion  rate  calculated  from  the 
reformate  gas  composition  exceeded  99.5%  at  reformer  outlet  tem¬ 
peratures  above  1013  K.  At  an  outlet  temperature  of  923 1<  a  fuel 
conversion  rate  of  86.2%  only  could  be  achieved. 

For  every  operation  point  of  the  catalyst  test  program  a  refor¬ 
mate  gas  sample  was  analysed  by  gas  chromatography  (GC)  after 
condensation  at  278  K.  The  GC  allows  the  detection  of  gaseous 
hydrocarbons  with  carbon  numbers  up  to  5.  In  Fig.  8  (left)  the 
detected  species  are  described  for  a  reformer  temperature  variation 
at  a  SCR  equal  to  4.  At  reformer  outlet  temperatures  above  1013  K 
only  methane  was  detected  in  the  reformate  gas.  The  concentration 
increases  with  decreasing  temperature  as  the  reforming  reaction 
is  constrained  and  the  fuel  conversion  rate  falls.  At  923  K  ethane, 
ethene,  and  propane  occurred  additionally  which  agrees  with  the 
observation  made  in  Fig.  7  that  total  hydrocarbon  concentration 
was  much  higher  or  conversion  rate  was  low  at  that  tempera¬ 
ture.  The  second  diagram  in  Fig.  8  (right)  shows  the  hydrocarbon 
species  detected  at  a  constant  set  temperature  equal  to  1013  K  and 


Fig.  8.  Left:  dry  reformate  gas  concentrations  at  SCR  4;  right:  dry  reformate  gas  concentrations  at  set  temperature  1013  K. 


Table  2 

Validation  cases,  conditions. 
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Case 

1  [16] 

2  [16] 

3  (Section  3.1) 

Pth>  reformer  inlet 

lOkW 

41<W 

61<W 

Pth,  burner  inlet 

4.25  kW 

2.5  kW 

3.8  kW 

SCR 

4 

4 

4.4 

Air  ratio,  k 

2.25 

2 

3.4 

Inlet  temperature  reformer 

873  K 

640 1< 

733  K 

Inlet  temperature  burner 

573  K 

573  K 

613  K 

Fuel  reformer 

Shellsol  surrogate 

Shellsol  surrogate 

Low  sulphur  fuel  oil,  6  wt.  ppm  sulphur,  long  term  test 

Catalyst  activity,  a 

100% 

100% 

3% 

a  variation  of  the  SCR.  Mainly  methane  was  detected  whose  con¬ 
centrations  increased  with  decreasing  SCR.  At  a  SCR  equal  to  3  also 
ethane  occurred  in  the  reformate  gas  with  a  concentration  equal  to 
32  vppm.  Based  on  these  results  the  lower  operational  limit  value 
for  the  SCR  can  be  set  to  3.5  and  for  the  reformer  temperature  to 
1013  K. 


5.  Computational  results  and  discussion 

5.1.  Validation 

Before  starting  with  the  computation  of  design  variations  CFD 
results  of  calculations  are  compared  to  experimental  results  of  the 
lOkW  reformer.  Fig.  9  shows  the  corresponding  outline  of  the  val¬ 
idation  case  with  the  position  of  inlets  and  outlets. 

In  Table  2  the  conditions  of  the  regarded  validation  cases  are 
listed.  Case  no.  1  and  2  were  measured  by  Wruck  [16].  Reformer 
reactant  was  Shellsol  surrogate  as  mentioned  in  Section  3.2.  Since 
the  fuel  contains  less  than  £s  <  2 wt. ppm  sulphur  catalyst  activity 
a,  given  by  the  following  equation,  remains  unaffected: 

a  =  |-  (12) 

->0 

Assuming  an  uniform  catalyst  activity  and  uniform  deactivation  the 
rate  constant  k0  related  to  the  mass  of  catalysts  are  corrected  with 
the  catalyst  activity  a  is  given  by: 

k*0  =  k0a  (13) 
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Fig.  9.  Outline  of  current  reactor. 


The  catalyst  activity  was  determined  indirectly  by  computation  of 
variations. 

In  all  cases  the  SCR  of  the  reformer  reactants  is  about  4  and  a  low 
sulphur  fuel  oil  containing  £s  =  5-6wt.ppm  sulphur  was  supplied 
to  the  burner.  In  contrast  to  the  first  two  variations  a  low  sulphur 
fuel  oil  was  used  in  the  reformer.  Beyond  that  the  air  ratio  was 
increased  to  A  =  3.4.  Due  to  the  fact  that  the  experiment  no.  3  (see 
Section  3.1)  was  a  long  term  test  with  low  sulphur  fuel  oil  and 
with  an  unadapted  regeneration  strategy  the  catalyst  activity  falls 
to  3%.  Beyond  that  the  used  catalyst  was  not  optimised  for  a  fuel 
containing  sulphur. 

Fig.  10  shows  the  temperature  distribution  in  the  plate  and  fuel 
mass  fractions  in  the  reformer  and  burner  for  case  no.  1.  Especially 
at  low  air  ratios  the  simple  lateral  inlet  for  the  burner  leads  to  a  hot 
outer  region.  Since  the  reaction  rate  of  the  steam  reforming  reaction 
is  supported  by  high  reactor  temperature  only  in  the  middle  of  the 
reactor  plate  hydrocarbon  reaches  the  outlet.  The  amount  of  fuel 
at  the  outlet  corresponds  to  a  fuel  conversion  of  X=92%.  Due  to  a 
relative  low  reactor  outlet  temperature  of  T=  973  K  a  full  conversion 
in  the  reformer  was  not  achieved. 

The  computational  results  are  validated  with  the  measured  tem¬ 
perature  distribution  in  the  reactor  and  fuel  conversion  in  the 
reformer  shown  in  Fig.  11.  In  all  cases  the  same  fuel  conversion 
was  reached.  The  agreement  of  the  temperatures  on  the  central 
line  is  very  good  for  case  1  and  3.  In  case  2  a  maximum  deviation 
of  23  K  can  be  observed.  Considering  the  error  of  thermocouples 
this  result  is  satisfactory.  The  temperatures  on  the  outer  lines  are 
also  in  a  good  agreement  with  a  maximum  absolute  deviation  of 
24  K.  First  calculation  with  a  more  complex  reaction  mechanism 
for  the  combustion  with  homogenous  reactions  a  temperature  dis¬ 
tribution  closer  to  the  measured  values  can  be  found.  The  measured 
values  at  the  burner  inlets  differ  in  some  case.  One  possible  reason 
is  a  different  coating  at  these  places.  On  grounds  of  low  catalyst 
activity  in  case  3  a  conversion  of  just  X=  67%  was  observed. 

5.2.  Variation  of  geometry 

The  aim  of  the  computational  investigation  is  the  optimisation 
of  performance  and  durability  of  the  reactor.  The  following  criteria 
are  stated: 

•  Fast  overheating  of  reformer  reactants  to  reduce  coke  formation. 

•  High  fuel  conversion  rates  in  the  reformer. 

•  Homogenous  temperature  distribution  to  prevent  material  dam¬ 
ages. 

•  Maximum  temperature  below  operating  temperature  limit 
T=  1173  K. 

•  Low  pressure  loss  in  the  burner  to  increase  system  efficiency. 

In  the  framework  of  the  numerical  investigation  the  following 
geometry  variations  are  rated  on  the  basis  of  the  computational 
results  of  case  3  and  case  1  with  increased  outlet  temperature 
T  =1023 1<: 
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Fig.  10.  Computation  results,  current  reactor,  case  1,  fuel  mass  fraction  reformer  (left)  and  burner  (middle),  temperature  in  the  plate  (right). 


Distance  from  reformer  inlet  /  mm  Distance  from  reformer  inlet  /  mm 


Fig.  11.  Comparison  of  temperature  distribution  in  reactor,  current  reactor,  cases  1-3,  measurements  and  computational  results,  central  line  (left)  and  outer  lines  (right). 


•  Current  reactor. 

•  Current  reactor  with  increased  channel  height. 

•  Current  reactor  with  increased  burner  inlet  width. 

•  Current  reactor  with  reduced  burner  inlet  width. 

•  Current  reactor  with  switched  sides  (reformer  and  burner). 

•  S-flow  reactor. 

•  Cross  flow  reactor. 

•  Counter  flow  reactor  based  on  the  current  reactor. 

•  Parallel  flow  reactor  with  supply  of  burner  reactants  through 
multiple  holes. 

•  Parallel  flow  reactor  with  switched  sides  (reformer  and  burner). 

Fig.  12  shows  the  computational  results  of  the  cross  flow  reac¬ 
tor.  Since  the  heat  demand  and  heat  source  do  not  match,  a  hot  spot 
appears  at  top  left.  The  standard  deviation  of  solid  phase  tempera¬ 
ture  for  case  3  and  the  cross  flow  reactor  is  100K  (see  Fig.  13).  The 
calculation  of  the  current  reactor  and  the  same  conditions  result 


in  a  deviation  of  just  60  K.  A  standard  deviation  of  81  K  and  117  K 
is  observed  for  S-flow  and  the  counter  flow  reactor.  Furthermore, 
operating  temperature  limit  is  exceeded  for  counter  flow  and  cross 
flow  reactor.  Table  3  presents  a  rating  matrix  regarding  the  geom¬ 
etry  variations.  The  current  reactor  with  switched  sides  has  a  good 
performance  with  conditions  referring  to  case  3,  but  at  low  air  ratios 
as  in  case  1  the  maximum  temperature  is  exceeded  and  standard 
deviation  of  91  K  is  reached.  In  the  current  reactor  and  for  case  3 
the  reformer  reactants  are  heated  up  to  T=  973 1<  in  the  range  up  to 
1 5  mm  of  the  reactor  length.  By  increasing  the  burner  inlet  width 
by  50%  the  length  for  reactants  overheating  rises  by  40%  and  vice 
versa.  So,  smaller  inlets  are  favoured  for  reactant  overheating.  By 
increasing  channel  height  from  0.6  mm  up  to  0.8  mm  the  pressure 
loss  in  the  burner  is  reduced  by  55%. 

For  high  air  ratios  referring  to  case  3  the  parallel  flow  reac¬ 
tor  with  supply  of  burner  reactants  through  multiple  holes  has  no 
advantages  over  the  current  reactor.  Rather  than  that,  the  length  for 
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Fig.  12.  Computational  results  cross  flow  reactor,  case  3,  outline  (left),  temperature  in  the  plate  (middle),  fuel  mass  fractions  reformer  (right). 
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Fig.  13.  Standard  deviation  of  temperature,  case  3,  geometry  variations. 


Table  3 

Rating  matrix  geometry  variations:  0,  does  not  meet  demands;  10,  meet  demands  perfectly;  0*,  exceeded  temperature  limit  lead  to  drop  of  rating;  C3,  case  3;  Cl,  case  1  with 
increased  reactor  outlet  temperature  1023 1<  and  reduced  catalyst  activity  3%. 


Criteria 

Overheating 

reformer 

reactants 

Conversion 

reformer 

Homogeneity 

temperature 

Maximum 

temperature 

Pressure  loss 
burner 

Total 

Weighting 

0.3 

0.35 

0.2 

0.1 

0.05 

1 

Conditions 

C3 

Cl 

C3 

Cl 

C3 

Cl 

C3 

Cl 

C3 

Cl 

C3 

Cl 

Current  reactor 

7 

9 

1 

4 

5 

4 

8 

1 

7 

9 

4.6 

5.5 

Increased  channel  height 

7 

9 

1 

4 

4 

5 

9 

2 

10 

10 

4.7 

5.8 

Increased  inlet  width 

6 

- 

1 

- 

5 

- 

8 

- 

8 

- 

4.3 

- 

Smaller  inlet  width 

8 

10 

1 

4 

4 

5 

8 

2 

5 

8 

4.6 

6 

Switched  sides 

7 

10 

2 

3 

6 

1 

8 

0 

8 

9 

5.2 

0* 

S-flow 

6 

- 

1 

- 

2 

- 

8 

- 

2 

- 

3.5 

- 

Cross  flow 

3 

- 

5 

- 

0 

- 

0 

- 

10 

- 

0* 

- 

Counter  flow 

5 

- 

9 

- 

0 

- 

0 

- 

6 

- 

0* 

- 

Parallel  flow 

5 

9 

1 

4 

3 

5 

8 

4 

8 

10 

3.7 

6 

Parallel  flow,  switched  sides 

5 

7 

1 

3 

2 

6 

8 

4 

6 

7 

3.4 

5.1 

reactant  overheating  increases  by  57%  as  shown  in  Fig.  14.  The  con¬ 
version  of  fuel  is  reduced  by  3%.  However,  at  low  air  ratio  (case  1) 
the  maximum  solid  phase  temperature  in  the  parallel  flow  reactor 
is  below  the  value  in  the  current  reactor.  Therefore  the  parallel  flow 
reactor  permits  lower  air  ratios  without  exceeding  the  operating 
temperature  limit. 

In  contrast  to  the  current  reactor  the  complexity  and  costs  of 
the  parallel  flow  reactor  is  higher. 


Fig.  14.  Calculated  temperature  distribution  in  the  reactor,  current  reactor  and  par¬ 
allel  flow  reactor,  case  1  and  case  3. 


Therefore  additional  CFD  simulations  are  performed  to  compare 
these  two  parallel  flow  concepts  under  different  conditions.  For 
this  a  variation  of  one  parameter  in  each  case  is  put  into  effect.  The 
conditions  at  reference  point  are  given  in  Table  4.  An  idealised  2D 
parallel  flow  reactor  was  used  for  this  computational  study. 

Fig.  15  shows  the  result  of  these  parameter  variations.  By 
increasing  the  outlet  temperature  via  thermal  input  in  the  burner 
the  level  of  temperature  is  increased  in  the  whole  reactor  for 
both  parallel  flow  concepts  (see  Fig.  15a).  However,  only  in  the 
current  reactor  the  operating  temperature  limit  is  exceeded.  For 
increased  air  ratio  A  =  3  the  temperature  distribution  in  the  reactor 
are  homogenised  in  both  cases  (see  Fig.  1 5b).  In  contrast  to  the  pre¬ 
vious  computations  the  length  of  superheating  reformer  reactants 
in  the  parallel  flow  reactor  is  shorter  than  in  the  current  reactor.  The 
difference  between  these  calculations  is  the  higher  catalyst  activity 
in  the  reformer  and  the  idealised  parallel  flow  reactor.  By  decreas¬ 
ing  the  catalyst  activity  the  reactor  temperatures  rise  about  a  fix 


Table  4 

Conditions  at  reference  point. 


Pth,  reformer  inlet 

lOkW 

Reactor  outlet  temperature 

1023  K 

SCR 

4 

Air  ratio,  X 

2 

Inlet  temperature  reformer 

773  K 

Inlet  temperature  burner 

573  K 

Catalyst  activity,  a 

100% 
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Distance  from  reformer  inlet  /  mm 


Distance  from  reformer  inlet  /  mm 


Fig.  15.  Parameter  variation  calculation,  current  reactor  and  ideal  parallel  flow  reactor,  variation  of  reactor  outlet  temperature  (top  left,  (a)),  variation  of  air  ratio  (top  right, 
(b)),  variation  of  catalyst  activity  (bottom  left,  (c)),  variation  of  thermal  power  supply  reformer  (bottom  right,  (d)). 


value  (see  Fig.  15c).  A  lower  reformer  inlet  power  leads  to  lower 
reactor  temperatures  (see  Fig.  15d).  At  the  reference  point  conver¬ 
sion  rates  of  Xr  =  100%  and  Xr  =  97%  are  calculated  for  the  parallel 
flow  reactor  and  the  current  reactor.  By  decreasing  the  outlet  tem¬ 
peratures  to  Tout  =  923%.  The  conversion  rate  drops  to  X!  =  93%  for 
the  current  reactor.  This  effect  was  also  observed  in  the  frame¬ 
work  of  the  experiments  performed  on  the  small  scale  reactor  (see 
Section  3.2). 

The  observed  hot  spots  in  the  current  concept  were  also  theo¬ 
retically  examined  by  Gritsch  [1].  He  investigated  the  conduction 
of  flow  in  catalytically  coated  plate  heat  exchanger  for  coupled 
methane  steam  reforming  and  methane  combustion.  During  com¬ 
putational  investigations  he  found  out  that  hot  spots  can  occur 
in  a  counter  flow  reactor.  In  a  parallel  flow  reactor  heat  source 
and  heat  sink  are  more  balanced.  Gritsch  mentioned  that  a  mul¬ 
tistage  supply  of  fuel  on  the  combustion  side  enhances  the  reactor 
performance. 

The  unfavourable  temperature  distribution  of  a  counter  flow 
reactor  is  confirmed  by  this  work.  Furthermore,  the  result  of  the 
numerical  study  is  that  two  different  parallel  flow  concepts  best 
meet  the  demands.  Referring  to  the  stated  requirements  for  the 
performance  of  a  steam  reformer  the  parallel  flow  reactor  with 
multiple  holes  and  the  current  reactor  provide  a  more  homogenous 
temperature  distribution  and  a  faster  reformer  reactants  over¬ 
heating  than  the  cross  flow,  S-flow  and  counter  flow  concepts. 
Beyond  that  the  temperature  distribution  in  the  parallel  flow  reac¬ 
tor  depends  on  the  following  conditions: 

•  Design  of  parallel  flow  reactor. 

•  Catalyst  activity. 

•  Thermal  input. 

•  Air  ratio. 

The  current  reactor  with  increased  channel  height  and  the  reactor 
with  reduced  burner  inlet  width  also  achieve  the  best  rating  results. 


The  reasons  are  the  reduced  pressure  loss  in  the  first  case  and  the 
enhanced  reactants  overheating  the  second  case. 


6.  Conclusions 

•  Experimental  tests  were  performed  on  a  lOkW  steam  reformer 
and  a  small  scale  reformer  for  fuel  oil  and  diesel  fuel. 

•  The  experimental  results  show  that  the  reformer  design  applied 
for  the  validation  allows  inlet  temperatures  below  773 1<  because 
of  its  internal  superheating  capability. 

•  CFD  computation  reproduces  the  measured  temperature  field  in 
the  reactor. 

•  Comparable  results  (temperature  profiles)  are  presented  for  the 
numerical  and  experimental  data  considering  three  different 
operating  points. 

•  The  aim  of  geometry  optimisation  is  enhanced  reactant  over¬ 
heating  of  reformer  reactants,  conversion  of  fuel  in  the  reformer, 
homogeneity  of  solid  phase  temperatures  and  lower  pressure  loss 
in  the  burner. 

•  The  present  reactor  (parallel  flow  with  lateral  inlets)  with 
reduced  burner  inlet  width  or/and  increased  channel  height  best 
meets  demands. 

•  For  low  air  ratios  the  parallel  flow  reactor  shows  an  improved 
performance  in  comparison  to  the  reference  system.  However, 
in  some  cases  the  time  for  reactant  overheating  on  the  reformer 
side  increases. 
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